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Recently Grassberger has shown that the correction to scal- 
ing for the conductance of a bond percolation network on a 
square lattice is a nonmonotonic function of the linear lattice 
dimension with a minimum at L = 10, while this anomalous 
behavior is not present in the site percolation networks. We 
perform a high precision numerical study of the bond percola- 
tion random resistor networks on the square, triangular and 
honeycomb lattices to further examine this result. We use 
the arithmetic, geometric and harmonic means to obtain the 
conductance and find that the qualitative behavior does not 
change: it is not related to the shape of the conductance dis- 
tribution for small system sizes. We show that the anomaly at 
small L is absent on the triangular and honeycomb networks. 
We suggest that the nonmonotonic behavior is an artifact of 
approximating the continuous system for which the theory 
is formulated by a discrete one which can be simulated on a 
computer. We show that by slightly changing the definition of 
the linear lattice size we can eliminate the minimum at small 
L without significantly affecting the large L limit. 

PACS numbers: 64.60.Fr, 64.60. Ak, 05.70.Jk 



I. INTRODUCTION 

According to finite size scaling theory, the conducance 
of a finite random resistor network is expected to vary 
with the system size: 

<a>=L-'[a + bf{L)] (1) 

Here t ~ 0.982 is the conductivity exponent and f{L) 
is the correction to scaling term which vanishes as the 
size of the system becomes infinite [^|-Q.In this work we 
study the behavior of the correction-to-scaling for finite 
systems that can be simulated on a computer. 

Recently Grassberger ||] numerically studied the con- 
ductance of random resistor networks on the square lat- 
tice by numerical simulations for 2 < L < 4096. He 
found that for bond percolation networks the corrections- 
to-scaling are nonmonotonic: there is a dip in the 
corrections-to-scaling present at small lattice sizes. Even 
though this effect is rather weak, it can be clearly distin- 
guished by high precision numerical simulations. In the 
site percolation networks this behavior is absent. Grass- 
berger showed that the bond percolation data could be 
fitted to the log-periodic form < cr > L* ~ sin(log(L)) 
and pointed out that only the comparison with site perco- 
lation resistor networks enables us to reject this. Because 
this behavior is not understood we wanted to investigate 



it further by studying systems defined on different lat- 
tices. One question we ask is whether this behavior is the 
result of the strong shape dependence of the probability 
distribution of conductances for small system sizes. To 
answer this we simulated square lattice percolating net- 
works and calculated the arithmetic, geometric and har- 
monic means (Grassberger |^ used the arithmetic mean 
to obtain the conductance) . Another question is whether 
the oscillations are universal or particular to the square 
lattice used in |]l|. To test the universality we also sim- 
ulated bond percolation resistor networks on triangular 
and honeycomb lattices. 

Our paper is organized as follows: in Section 2 we give 
the details of our numerical work and present our results. 
Section 3 contains the discussion and conclusions. 

II. NUMERICAL SIMULATIONS 

The systems we study are random resistor networks 
that are connected to perfectly conducting busbars at 
the two opposite edges of the array, while the transverse 
boundaries are free. To simulate a random resistor net- 
work we start with a bond percolation network and assign 
unit conductance to the bonds that are present while a 
missing bond gives zero conductance. The calculation of 
the network conductance for a given configuration of con- 
ductances is done by the Lobb- Frank algorithm Q , which 
uses a succession of star-triangle (ST) and triangle-star 
(TS) transformations to reduce the network to an equiv- 
alent chain of resistors. 

Our square lattice simulations were performed on rect- 
angles {0 < X < L,l < y < L), with the perfectly 
conducting busbars placed at a; = and x = L. This 
means that the lattice (busbars excluded) has dimen- 
sions {L — 1,L). The reason for this choice is to en- 
sure that the bond system is self-dual (see Ref. 3 and 
references therein). To simulate systems on a triangu- 
lar lattice we add a diagonal resistor to each unit cell 
of the square lattice without changing the position of 
the busbars. This way of representing the triangular net- 
work changes the geometry of the system - effectively the 
shape of the boundary becomes a rhombus |^ . Simula- 
tions on the honeycomb lattice arc done by starting with 
a honeycomb resistor network and replacing the leftmost 
column of resistors with perfect conductors and attach- 
ing the rightmost column of nodes to a perfectly con- 
ducting busbar. The pictures of these networks can be 
found in Ref. 4. The systems were taken to be at the 
percolation threshold, where Pc = 0.5 for square lattice, 
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Pc = 2sin(7r/18) = 0.34729 for the triangular lattice and 
Pc = 1 ~ 2sin(7r/18) = 0.65271 for the honeycomb lat- 
tice For the square lattice the conductance was 
obtained by exact enumeration for L < 4, while for the 
triangular lattice exact enumeration was used to obtain 
the result for L = 2. For other system sizes the number 
of random configurations generated was at least 10''. 

Our results for the square lattice are shown in Fig. 1, 
for the triangular lattice in Fig. 2 and for the honeycomb 
lattice in Fig. 3. We plot the conductance multiplied by 
L raised to the conductivity exponent t = 0.982 as a 
function of L on a semilogarithmic scale. 
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FIG. 1. Average conductance on square lattices of size 
Lx L, multiplied by L"'®*^. The upper set of data is obtained 
by using the arithmetic mean, the middle set is obtained from 
the geometric mean and the lower set is the harmonic mean. 
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FIG. 2. Average conductance on triangular lattices of size 
L X L, multiplied by L"'^^^ . The symbols have the same 
meaning as in Fig. 1 
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FIG. 3. Average conductance on honeycomb lattices of size 
L X L, multiplied by L"'^*^. The notation is the same as in 
Fig. 1 and Fig. 2. 

In these figures full points (circles) are obtained by 
using the arithmetic mean, stars are obtained from the 
geometric mean and crosses correspond to the harmonic 
mean (where the averages are computed using only the 
samples with nonzero conductance). The errors are 
smaller than the size of the points on the graphs. The 
points in Fig. 1 obtained by the arithmetic mean match 
the data obtained by Grassberger in |^]. 



III. DISCUSSION AND CONCLUSIONS 

Let us first consider our results for the networks defined 
on a square lattice. 

The data obtained by using the arithmetic mean start 
at a low value for the lowest lattice size {L = 2) and 
then there is a jump to a high value at lattice size L = 3 
followed by the monotonic decrease up to L « 10. If 
we use the harmonic or geometric mean the jump is not 
present - < a > L* is a monotonically decreasing function 
until it reaches a minimum. Thus the behavior at L = 2 
is affected by the peculiar (and very discrete) distribution 
of conductances. 

In all three cases a minimum is present at low L. Our 
numerical investigation of conduction on a square lattice 
shows that the dip in the correction-to-scaling term for 
small systems does not get averaged out by using different 
means to obtain the conductance. We did not find the 
anomalous behavior for the systems defined on triangular 
and honeycomb lattices. The dip at small L is not present 
and the behavior is quite similar to those obtained in Ref. 
1 for the case of site percolation on a square lattice. 

The anomaly present in the case of square lattice is 
rather weak - only high precision simulations are able to 
reveal it. We notice that changing i by a fraction of the 
lattice spacing will strongly affect the behavior of < a > 
i* for small L while the large L results will be practically 
unaffected. Replacing L by L — A as the independent 
variable gives a monotonic function for A = 0.05; for 
A = 0.1 the corrections-to-scaling behave qualitatively 
in the same way as the results obtained for the other 
systems mentioned above (Ziff used a similar constant 
in connection with crossing probability problems). This 
is illustrated in Fig. 4. 
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FIG. 4. Average conductance on square lattices of size 
L X L, multiplied by (L - A)"'^®^. The upper set of data is 
obtained by using the arithmetic mean (this set is represented 
in Fig. 1. by the same symbol), the middle set is obtained 
by correcting L by 0.05 and the lower set is obtained with a 
correction of 0.1. 

In the case of triangular and honeycomb systems the 
rescaling of L does not qualitatively change the behavior 
of < (T > L* - in particular, the oscillation at small L 
does not occur. This is shown in Fig. 5 and Fig. 6 for 
triangular and honeycomb lattices, respectively. 
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FIG. 5. Average conductance on triangular lattices of size 
L X L, multiplied by (L - A)°-^*^. The middle set of data is 
obtained by using the arithmetic mean (this set is represented 
in Fig. 2. by the same symbol), the lower set is obtained by 
correcting L by 0.05 and the upper set is obtained with a 
correction of -0.05. 
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FIG. 6. Average conductance on honeycomb lattices of size 
Lx L, multiplied by (L - A)°-^®^. The middle set of data is 
obtained by using the arithmetic moan (this set is represented 
in Fig. 3. by the same symbol), the lower set is obtained by 
correcting L by 0.05 and the upper set is obtained with a 
correction of -0.05. 



The justification for this adjustment of L is that the 
theory is formulated in the continuum hmit while the dis- 
crete systems we simulated have the smallest scale equal 
to the lattice spacing. While for the large systems this 
should not matter, the small system results might be af- 
fected by this if we are doing very precise simulations. 

The author thanks J. P. Straley for critical reading of 
the manuscript. 
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